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Abstract 

In this paper, a systematic analysis of hourly wind speed data obtained from four potential wind generation 
sites in North Dakota is conducted. The power spectra of the data exhibited a power law decay characteristic 
of 1/ f a processes with possible long range correlations. The temporal scaling properties of the records were 
studied using multifractal detrended fluctuation analysis MFDFA. It is seen that the records at all four 
locations exhibit similar scaling behavior which is also reflected in the multifractal spectrum determined 
under the assumption of a binomial multiplicative cascade model. 

1 Introduction 

Wind energy is a ubiquitous resource and is a promising alternative to meet the increased demand for energy 
in recent years. Unlike traditional power plants, wind generated power is subject to fluctuations due to the 
intermittent nature of wind. The irregular waxing and waning of wind can led to significant mechanical stress 
on the gear boxes and result in substantial voltage swings at the terminals, £Q. Therefore, it is important to 
build suitable mathematical techniques to understand the temporal behavior and dynamics of wind speed for 
the purposes of modeling, prediction, simulation and design. Attempts to identify the features of wind speed 
time series data were described in |2] and [3]. To our knowledge, the first paper to bring out an important 
feature of wind speed time series was 0] . In ^ , the authors examined long term records of hourly wind speeds 
in Ireland and pointed out that they exhibited what is known as long term dependence. Seasonal effects, spa- 
tial correlations and temporal dependencies were incorporated to build suitable estimators. Evidence for the 
presence of long memory correlations was provided by inspecting the periodogram of the residuals from a fitted 
and autoregressive model of order nine i.e. AR(9) 0]. 

For the pursuit of wind energy development in the state of North Dakota which has an estimated potential of 
250 GW, j^j, several wind monitoring stations have been set up where the hourly average wind speeds along 
with meteorological variables are being recorded. As the power developed by the wind turbines is proportional 
to the cube of the wind speed, the recorded data can be used to assess the wind power potential of the site. In 
addition, an analysis of the recorded data can be useful in understanding the nature of the fluctuations in power 
developed by the wind turbine. Conventionally, the statistical modeling of wind speed data has been done by 
using Weibull and Rayleigh probability distribution functions, [Sj- While these methods are useful in providing 
estimates of wind generated energy yield, they do not explicitly bring out the nature of the temporal variations in 
wind speed and thus, wind power. Complex motions in the atmosphere tend to render the wind speed distinctly 
non-stationary and intermittent. Traditional power spectral (or equivalently the auto correlation function) 
analyses are restricted in their scope for the present study due to their susceptibility to non-stationarity. Power 
spectral techniques have been useful successfully to detect possible long-range correlations of the form of the 
form S(f) ~ l//' 3 , [7j. Long range correlations generally indicate that temporally well separated samples of the 
time series are correlated with each other and indicative of self-similar (sclf-affine) behavior. Self similar time 
series can be characterized by 
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y(t) ee a a y(t/a) (1) 

where the = in Eqn.Q denotes that both sides of the equation have identical statistical properties. The 
exponent a in Eqn.JQ) is called the self-similarity parameter, or the scaling exponent. Values of a in the 
range (0, 0.5) characterize anti-persistence, whereas those in the range (0.5,1) characterize persistent long range 
correlations with a = 0.5 representing uncorrelated noise. The exponent (/?) estimated from the power spectrum 
(S(f) ~ 1/ f 13 ) is related to a as /3 = 2a — 1. The temporal trace and corresponding power spectrum of one of 
the representative records is shown in Fig. ^ 




Figure 1: Temporal trace and Power spectrum of a representative record 

Examining the power spectrum of the wind speed indicates a power-law decay of the form S(f) ~ l//' 3 . The 
auto-correlation functions (ACF) decay slowly to zero and the first zero crossing of the ACFs occur at lags of 
61,56,60 and 57 respectively for the four data sets. Such features are characteristic of statistically self similar 
processes with well defined power law correlations, [7]. 

Estimating the Hurst exponent (H) from the given data is an alternate effective way to determine the nature 
of the correlations in it [8] . Hurst exponents have been successfully used to quantify long range correlations in 
plasma turbulence [5], ^U], finance network traffic, ^3] an d physiology, Long range correlations 

are said to exist if 0.5 < H < 1 (see ^5], ^B] for details). There are several methods such as the variance 
method, re-scaled range (R/S) analysis, Whittle estimator and wavelet based methods to estimate the Hurst 
exponent. Hurst estimators arc susceptible to such artifacts such as polynomial trends which cannot be ruled 
out in experimental data and therefore, they may give spurious results. Recently the Detrended Fluctuation 
Analysis (DFA), ^7J and its extensions have been proposed as an alternate effective technique to determine 
possible long range correlations in data sets obtained from a diverse settings, . |20|. |21j . 

It should be noted that the techniques listed above can only extract a single scaling exponent from a time series. 
However, it is possible that the given process may be governed by more than one scaling exponents, in which 
case a single scaling exponent would be unable to capture the complex dynamics inherent in the data. Therefore, 
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these methods are appropriate only for the analysis of monofractal signals which have uniform scaling properties 
throughout the signal which can be characterized by a single exponent. On the other hand, multifractal signals 
are far more complex than monofractal signals and require more than one (theoretically infinite) exponent to 
characterize their scaling properties, |14|. 

In the present study, the wind speed time series data is studied using a fairly robust and powerful technique 
called Multifractal Detrended Fluctuation Analysis (MFDFA), |22| . The method provides a systematic means to 
identify and more importantly quantify the multiple scaling exponents in the data |22| . The scaling exponents 
of the data are estimated under the assumption of a binomial multiplicative cascade model. This is carried out 
on the wind speed data acquired from four spatially separated monitoring locations in North Dakota. The paper 
is organized as follows. In Sec[3 the acquisition of wind speed data is described. In Sec. |3| a brief description 
of multifractal detrended fluctuation analysis is provided. In SecQ] the results of MFDFA and the calculation 
of the multifractal spectrum is provided. The conclusions are summarized in SecJS] 

2 Data Acquisition 

The hourly averaged wind speeds at four geographically well separated locations in North Dakota were recorded 
by means of a conventional cup type anemometers located at a height of 20 m. The co-ordinates of the map 
are provided in Tabled From Fig[3 it can be seen that the locations are sptially well separated across North 
Dakota. The wind speed data to be discussed were recorded over a time frame ranging from 11/29/2001 to 
07/28/2003. 



Table 1: Co-ordinates of the wind monitoring sites 



Station 


Latitude 


Longitude 


Elvn (ft) 


Site 1 


N 47 27.84' 


W 99 8.18' 


1570 


Site 2 


N 46 13.03' 


W 97 15.10' 


1070 


Site 3 


N 48 52.75 


W 103 28.4' 


2270 


Site 4 


N 46 12.67' 


W 103 13.07' 


2880 




Figure 2: A map of the wind monitoring sites 



3 Multifractal Detrended Fluctuation Analysis (MFDFA) 

MFDFA, |221 a generalization of DFA has been shown to reliably extract more than one scaling exponent from 
a time series. A brief description of the algorithm is provided in Section III, for completeness. A detailed 
explanation can be found elsewhere |22]. Consider a time series {xk},k = 1...N. The MFDFA algorithm 
consists of the following steps. 
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1. The series {xk} is integrated to form the integrated series (also called "profile") [yu] given by 

i—k 

y(k) = 5>(i) - z] k = l,...N (2) 

where x represents the average value. 

2. The series is divided in to n s non-overlapping boxes of equal lengths where n s = int(N/s). To 
accommodate the fact that some of the data points may be left out, the procedure is repeated from the 
other end of the data set |22j . 

3. The local polynomial trend y v with order v is fit to the data in each box, the corresponding variance is 
given by 

i—s 

i=i 

for v = 1, .. .n s . Polynomial detrending of order m is capable of eliminating trends up to order m-1. |23| 

4. The qth order fluctuation function is calculated from averaging over all segments. 

{i=2n s \ 1 / q 

^- y:\f\v, S )r ,2 \ (4) 

In general, the index q can take any real value except zero. 

5. Step 3 is repeated over various time scales s. The scaling of the fluctuation functions F q (s) versus the 
time scale s is revealed by the log-log plot. 

6. The scaling behavior of the fluctuation functions is determined by analyzing the log-log plots F q (s) versus 
s for each q. If the original series {xk} is power-law correlated, the fluctuation function will vary as 

(5) 

3.1 Multifractal Analysis 

In the MFDFA procedure (Section III), the exponent h(q) describes the scaling behavior of the qth order fluctu- 
ation function. For positive values of q, h(q) describes the scaling behavior of segments with large fluctuations 
while those of negative values of q, describe scaling behavior of segments with small fluctuations, (221 - F° r sta- 
tionary time series, the exponent h(2) is identical to the Hurst exponent. Thus the exponent h(q) is called as the 
generalized Hurst exponent, |22j . For monofractal time series which are characterized by a single exponent over 
all scales, h(q) is independent of q, whereas for a multifractal time series, h(q) varies with q. This dependence 
is considered to be a characteristic property of multifractal processes, |22j . 

A multifractal description can also be obtained from considering partition functions ( see [21] for details) 

v=n s 

z q (s)=J2 l^-y ( ,-i) S | 9 ~s r(9) (6) 

where r(q) is the Renyi exponent. A linear scaling of r(q) with q is characteristic of a monofractal data, whereas 
a nonlinear scaling is indicative of multifractal behavior. The h(q) obtained from MFDFA is related to the the 
Renyi exponent r(q) by 

qh(q) = r(q) + 1 (7) 
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Therefore, another way to characterize a multifractal series is the multifractal spectrum f(oth) defined by, 7 

In (9), ah is the Holder exponent which characterizes the singularities (cusps, ridges, chirps, spikes) in a process 
X(t) at time t, The multifractal spectrum f(cth) describes the singularity content of the process, i.e. the 
dimension of the set of times t where = a, 25 . Substituting for r(q) from Eqn.Q, we get, 

a h = h(q) + q—^i f( a h) = q[ah - h(q)] + 1 (9) 
The generalized exponents h{q) in Eqn.@ can be estimated by the formula 

assuming a binomial multiplicative cascade model |7] which has served as one of the standard paradigms to 
describe multifractal scaling, (see 0, |22 f° r details). Finally, the width of the multifractal spectrum 
f(a h ) at / = 0, given by 

,/ \ , / \ (Mb) — ln(a)) ,„„, 
Aa h = ft(-oo) - h(oo) = y , o 11 

ln2 

is computed to compare the strength of multifractality for different records. Incidentally, the multifractal 
spectrum can also be estimated using wavelets by the wavelet transform modulus maxima (WTMM) method, 
|27| which we do not pursue here. 



4 Results 

The log- log plot of the fluctuation functions F q (s) vs s for the wind speed records at all four locations with (q = 
-6, -4, -2, 2, 4, 6) using fourth order polynomial detrending is shown in FigEJ The generalized hurst exponent 
h(q) estimated for (q = -6, -4, -2, 2, 4, 6) using Eqn.JSJ) is shown in FigQ] One can note that at all locations, 
the slopes h(q) decrease as the moment is increased from negative to positive values. For example, the slope 
decreases from 0.88 when q = —6 to 0.6989 when q = +6 for site 3 which is indicative of multifractal behavior. 
Recent studies have suggested the choice of constrained randomized shuffled surrogates to reject the null that 
the observed scaling is due to the distribution as opposed to the temporal structure in the data. Constrained 
randomized shuffles (surrogates) are similar to traditional bootstrapping without replacement, i.e. each member 
in the original data is retained in the surrogate. While the temporal structure is destroyed the distribution of the 
original data is retained in the surrogate realizations. The surrogate realizations are equivalent to uncorrelated 
noise. A good exposition of the concepts of surrogate analysis can be found in |28|. In the present study, the 
MFDFA studies are performed on random shuffle surrogates of the original data, to reject the claim that the 
observed scaling exponent obtained on the original data is due to its distribution and not due to correlations 
present in it. The results shown in Fig. [^indicate that for all values of q, the fluctuations scale as F q (s) ~ s ' 5 
which signifies a loss of correlations. In other words, there is a loss of multifractality because the slopes h(q) 
are independent of q. The parameters a and b determined under the assumption of a binomial multiplicative 
process Ean. (|10J) is also shown FigQ] Using the values of a and b with Eqns.JSJ), © and the multifractal 
spectrum at all four locations shown in Figgis computed. The width of the multifractal spectrum is calculated 
from Ean. lfTTJl . 



5 Conclusions 

In this study, the hourly wind speed records at four geographically well separated locations in North Dakota 
were examined. The focus was on determining the nature of the correlations and hence, the dynamics in the 
wind speed data. The power-law decay of the power spectrum suggested the presence of possible long range 
correlations and self-similar behavior. A number of techniques proposed in the past to quantify self-similarity 
in a given data have been susceptible to non-stationarity which can manifest itself as polynomial trends. This 
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Figure 4: Generalized scaling exponents for wind speed records at the four sites. Circles indicate the actual 
values of the scaling indices for q — —6, —4, —2, 2, 4 and 6. The solid line represents the fit from the multiplicative 
cascade model (Ean. ljlOp ) obtained from a nonlinear least squares procedure. The values of a and b are indicated 
for each site. 
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Figure 5: Fluctuation functions at all sites using MFDFA for random shuffle surrogates of the original data. 
Symbols used to indicate the various moments are q — — 6 (-), q = — 4(inverted triangle), q = — 2 (x), q = 2 
(.), q = 4 (o), q = 6 (+). The curves are vertically shifted for clarity. Note that the curves are parallel to each 
other with a slope of ~ 0.5. 




Figure 6: Multifractal spectra of wind speed records at the four sites 
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motivated us to use a fairly robust technique such as the MFDFA. Analysis using MFDFA elucidated the depen- 
dence of h(q) on q, which is a hall mark of multifractal processes. The binomial multiplicative cascade model 
has been widely used to determine the multifractal spectra. The multifractal spectra (/(a)) and the widths 
(Aq) (assuming a binomial multiplicative cascade) were estimated at the four locations. 

Temporal and spatial variations of wind speed are influenced by several factors such as pressure gradient, 
turbulence, temperature and topography. Our preliminary studies suggest that long term wind speed variations 
possess a multifractal structure. Moreover, despite the heterogeneities across the spatially separated locations, 
certain quantitative features such as the spectrum (/(a)) and multifractal widths (Aa) seem to be retained. A 
deeper investigation of these issues with extensive data sets over wider geographical regions may provide clues 
towards understanding the nature of long term atmospheric wind speed variations. 
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